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A method is developed to determine the eigenvalues and eigenfunction of two-boson 2x2 matrix 
Hamiltonians include a wide class of quantum optical models. The quantum Hamiltonians have 
been transformed in the form of the one variable differential equation and the conditions for its 
solvability have been discussed. We present two different transformation procedure and we show 
our approach unify various approaches based on Lie algebraic technique. As an application, solutions 
of the modified Jaynes-Cummings and two-level Jahn Teller Hamiltonians are studied. 
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INTRODUCTION 

The algebraic techniques have been proven to be useful in the description of the physical problems in a variety 
of fields. Recently a new algebraic approach, essentially improving both analytical and numerical solution of the 
problems, has been suggested and developed for some nonlinear quantum optical systems 0, 0, 0, Q- Most of such 
developments are mainly based on linear Lie algebras, but it is evident that there is no physical reason for symmetries 
■ to be only linear. Nonlinear Lie algebra techniques and their relations to the nonlinear quantum optical systems 
have been discussed 0, 0, S 0- I n both cases finite part of spectrum of the corresponding Hamiltonian can be 
exactly obtained in closed forms and these systems are known as quasi-exactly-solvable(QES) termed by Turbiner 
and Ushveridzejlfj . It has been proven that the single boson Hamiltonians also lead to a QES under some certain 
O ' constraints [Til Il2| . 

In recent years there has been a great deal of interest in quantum optical models which reveal new physical 
phenomena described by the Hamiltonians expressed as nonlinear functions of Lie algebra generators or boson and/or 
fermion operators |l3l Il4l lid lid Il7j. Such systems have often been analyzed by using numerical methods, because 
the implementation of the Lie algebraic techniques to solve those problems is not very efficient and most of the other 
analytical techniques do not yield simple analytical expressions. They require tedious calculations [Til H^. Eof . 

The aim of this paper is to determine solvability of the two-boson 2x2 matrix Hamiltonians and discuss their 
applications and possible symmetry groups. The results of our procedure include the solutions of the Hamiltonians 
that possessing su(2), su(l, 1), Sp(4,R), osp(2, 1) and osp(2, 2) symmetries. The procedure presented in this paper 
also leads to the constructions of the nonlinear Lie algebras |21|. The Hamiltonians discussed here are not only 
mathematically interesting but they have potential interest in physics |22t l23l l24f. In order to keep our discussion 
simple we concentrate our attention to the solution of the Hamiltonians that include two boson and one fermion. 
Therefore we provide a first step toward the extension of the technique to the solution of the multi boson fermion 
C3 ' systems. 

The paper is organized as follows: In section 2, we construct a Hamiltonian including two boson operators in 
an arbitrary order and one fermion operator. Solutions of the Hamiltonians by using the invariance of the number 
operator have been discussed. As a practical example, solution of the Modified Jaynes-Cummings and two level Jahn 
Teller Hamiltonians have been obtained. In section 3 we present two transformation procedure that is appropriate to 
determine the conditions of the (quasi)exact solvability of the Hamiltonian. In section 4 we discuss the symmetry prop- 
erties of the Hamiltonian. The results of our discussion are that the Hamiltonian contains various Lie (super)algebras, 
namely, su(2), su(l, 1), Sp(4, R), osp(2, 1) and osp(2,2). We also point out that nonlinear Lie algebras can be con- 
structed as spectrum generating algebras. As an application we visit Jaynes-Cummings Hamiltonian including Kerr 
nonlinearity. The importance of our approach is that it provides a unification of the various approaches. Finally, in 
section 5 we comment on the validity of our method and suggest the possible extensions of the problem. 
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TWO BOSON ONE FERMION HAMILTONIAN AND ITS DIFFERENTIAL REALIZATION 



Two mode bosonic 2x2 matrix Hamiltonians play an important role in nonlinear quantum optical systems. The 
Hamiltonians of such systems can be generalized as follows: 
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f+ Twi J 7ri2)W3 ) m4-^ ~t~ ^ ^ ^711,712,713.714-^ • (1) 

Tlli 71^ 

The bosonic operator i?^ reads as follows: 

b u = (4) ui (air 2 (4r(a2r (2) 

The constants on, (3i, ji and <5i are related to the parameters of the physical Hamiltonian and ki,£i, rm and determine 
order of the interaction. The boson creation, a\, a 2 , and annihilation a\, a 2 operators obey the usual commutation 
relations 

[oi,Oi] = [af, af] =0, [^4] = {J;.JJ (3) 



and er± are Pauli matrices: 



o\ /o n / 1 , 

1 ' a+ = ' a ° = -1 ) ' (4) 



The number operator of such system can be expressed in the form 

N = sa^ai + pa 2 a 2 + ro$ (5) 

and satisfy the commutation relations 

[N,af] = sa^, [N,ai] = —sai, [N.a^j—pa^, [N, a 2 ] = -pa 2 

[N,a ± ] = ±2ra± [N,a ]=0. (6) 
The action of the number operator on the state \n\, n 2 ) is given by 

N\ni,n 2 ) = (sm +pn 2 + r) \ni,n 2 ) (7) 
The eigenvalue equation © leads to the following solution: 

Im.na) = (afy<j>i (4(4)" 5 ) It) + K)^^^ (4(«iT f ) ID • (8) 
where ||) and |J,) are up and down states and j is given by 

P 

j = n 1 + -n 2 . (9) 

s 

If N and H commute, the function (JSJ is also eigenfunction of the H . Therefore it is worth to seek the conditions for 
the commutation of H and N. This can easily be done by using the commutation relations © and © and we obtain 
the following relation 

[N, H] = - fc 2 ) + p(fc 3 - k 4 )]a klk2k3ki B k + 

ki 

CT + X^I S ( mi ~~ " ? ' 2 ) +^( m 3 ~ m 4) + r ]7mim 2 m3m 4 B m + (10) 

o-- 2j[s(rii - n 2 ) +p(n 3 - n 4 ) - r]<5 nin2n3 „ 4 B™ 
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Then the constant of motion N and H commute when the fallowing set of equation is satisfied 

s{k\ - k 2 ) +p(h - fc 4 ) = 

s(£i-£2)+p(£s-h) = 
s(mi — m 2 ) +p(m3 — 1714) + 2r = (11) 
s{n\ — n 2 ) +p(ri3 — 714) — 2r = 0. 

Now we demonstrate application of the procedure on a physical example. We study modified Jaynes-Cummings 
Hamiltonian and this give us an opportunity to test our approach because those Hamiltonians have been studied in 
literature. 

The modified Jaynes-Cummings Hamiltonian 

The modified Jaynes-Cummings Hamiltonian have been constructed to investigate sin gle two level atom placed in 
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the common domain of two cavities interacting with two quantized modes. It is given by 

H = u(afax + a 2 a 2 ) + -y-"o + Ai(aier + + a^cr-) + A 2 (a 2 cr + + a 2 f(T -)- (12) 

When the parameters 

a i,jM.l = 0i,j,k,l = li,j,k,l = = (13) 

except that 

a u ° 
"1,1,0,0 — ao,o,i,i = W)Po,o,o,o — 

7o,i,o,o = <$i, 0,0,0 = Ai, 70,0,0,1 = <$o, 0,1,0 = A 2 (14) 

then the Hamiltonians and l|12l) are identical. The condition (| 1 1 1> is satisfied when s = p = 2r. Before going further 
in the following we use the Bargmann-Fock representation, where creation and annihilation operators are replaced by 
multiplication and differentiation operators: 

aj = Zi, <n = ^7 (15) 
with respect to complex variable Zi . The eigenfunction (jSJ) is of the form 

\ ni ,n 2 ) = ( Zl ycP(x) |T) + ( Zl y +1 cp{x) ||) (16) 

where x — {z\)~ x z 2 . The solution of this system describes a quantum mechanical state of H provided that <j>(x) 
belong to the Bargmann-Fock space. The scalar product should be complete and normalizable, 

J(x)(t){x)^ w ^ dx d{Rex)d(Imx) < 00 (17) 

where W{x) is the weight function. The eigenvalue equation of the modified Jaynes-Cummings Hamiltonian can be 
written as 

H\m,n 2 ) = E\n u n 2 ) (18) 
Insertion of l|15fl and (|16f) into (|18fl yield the following two set of differential equation 
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ju + — — E <pi{x) + (j + l)\ 2 <t> 2 (x) + (A 2 - x\x) 



d(t> 2 (x) 



dx 

{j + 1)0,-^-^1 02 (a;) + (Ax + x\ 2 )Mx) = i] 1 L!)| 
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Bargmann-Fock space solution of the (|19f) can easily be obtained and they are given by 

M*) = Co(X2-x\ 1 y- n+1 (x 1 + x\ 2 ) n - 1 

where n is an integer and eigenvalues of the Hamiltonian is given by 



(x) = C 1 (X 2 -xX i y- n+1 (X 1 +xX 2 ) n (20) 



E = \ ((2j + IV ± \pHX\ + X\) + (ljo ~ ^) 2 ) • (21) 

Consequently we have obtained exact result for the eigenvalues of the modified Jaynes-Cummings Hamiltonian. The 
same result have been obtained in|5^], in the framework of the su(2) algebra. The procedure given here can be applied 
to obtain eigenfunction and eigenvalues of the various physical Hamiltonians. In the following example we consider 
the solution of the Jahn- Teller distortion problem. 

Two-Level Jahn- Teller Distortion Problem 

The well-known form of the JT Hamiltonian describing a two- level fermionic subsystem coupled to two boson modes 
has been given by |2fij : 

H = a+ai + a+a 2 + 1 + (- + 2^)cr + 2k[(oi + a+)a + + (a+ + a 2 )<7_]. (22) 

Our task is now to demonstrate the Hamiltonian l|22|l can be solved in the framework of the procedure given pre- 
vious section. It will be shown that our approach relatively very simple when compared previous approaches. The 
Hamiltonians (1221 and Q are identical, when the parameters are constrained to: 

Oil, 1,0,0 = £*o, 0,1,1 — a 0, 0,0,0 = 1) A), 0,0,0 = (r + 2/i)> 

70,1,0,0 = £1,0,0,0 = 70,0,1,0 = #0,0,0,1 = 2k (23) 

otherwise 

a i,j,kJ = 0i,j,k,l = li,j,k,l = = (24) 

The condition is satisfied when s = —p = 2r. In the Bargmann-Fock space the eigenfunction JHJ takes the form 

\n u n 2 ) = { Zi y<h (x) |T) + (ziy-'fa (x) |I) . (25) 

where x = z\z 2 and j : = m — n 2 . Since number operator N and Hamiltonian i!22[) commute, they have the same 
eigenfunction. Substituting into if^l) we obtain the following set of equation: 

x^[(x)+KX^' 2 (x) + -(3-2E + 2j + 4fj l )<j> 1 (x)+K(l+j + x)^2(x) = (26) 

(x) + x<j)' 2 (x) + | (3 - 2E + 2j - 4n)(t> 2 (x) + nfa (x) = (27) 

These coupled differential equations represent the Schrodinger equation of the E g> e Jahn- Teller system in the 
Bargmann's Hilbert space. The equation is quasi-exactly-solvable and success of our analysis leads to the solu- 
tion of the various quantum optical systems. The physical systems described by the differential equations (|26II27|) are 
discussed in|27|. 

The validity of the procedure depends on the choice of the 0!j, A, 7j and <5i.One can easily be obtain various physical 
Hamiltonians by appropriate choice of a,, (3i, 7^ and Si and by considering the conditions given in (jlip. 
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TRANSFORMATION OF THE OPERATORS 



In this section we discuss transformation of the Pauli matrices and boson operators. These transformations plays a 
key role to construct (quasi)exactly solvable 2x2 matrix Hamiltonians. The transformation can be done by introducing 
the following similarity transformation induced by the metric 

S = ( a +y4°-i+<i<y+°- (28) 

where c and d are constants. Since a±, and a 2 commute and u± t o also commute with the bosonic operators, the 
transformation of a\ and a+ under S can be obtained by writing a+ = e b , with [a\, b] = [a+ , b] = 0, 

ScuS- 1 = ai(a+)~ c 

Sa+S- 1 = a+{a+f (29) 

the transformation of a 2 and is as follows 

Sa 2 S -1 = a 2 — (ca+ai + d<r + (j-)(a+)- 1 

Sa+S' 1 = a+ (30) 



and the transformations of the a± are given by 

Sa±S- 1 =<J ± (a + ) ±d . (31) 

Before constructing one variable (quasi) exactly solvable differential equation of the Hamiltonian 0) under the 
transformations of the bosonic and fermionic operators by S, let us consider the other transformation operator: 

T = ( a2 y4"i+n<r+<r- (32) 

where e and r\ are constants. By using the similar arguments given in the previous transformation operations one can 
easily obtain the following transformations: 

T ai T- x = ai{a + )' e 
Ta+T- 1 = a+(a 2 y 

T^T- 1 = a 2 (33) 
Ta+T -1 = a+ + (ea+ai + ija + a-)(a2)- 1 
Ta ± T- x = a ± (a+) ±Tl . 

These transformations leads to the construction of the various differential realizations of (JIJ , depending on the choice 
of c, n, e and r\. In the following we discuss the possible forms of the Hamiltonian Q and their solutions. 



S-TRANSFORMED HAMILTONIAN 



In this section we discuss the solvability of the Hamiltonian QJ. The number operator TV describe the states of the 
corresponding Hamiltonian. The transformation of N under the operator S is given by 

TV' = SNS -1 = (s~ pc)a+a! + pa+a 2 + (2r - pd)a + a^ - r (34) 

The Hamiltonian Q is characterized by the total number of a± and a 2 bosons when the number operator TV commute 
with the whole Hamiltonian. In the transformed case it is only the number of a 2 bosons that characterize the system 
under the condition c — s/p and d — 2r/p. When the representation is characterized by a fixed number a+a 2 = j, 
then the transformed form the Hamiltonian can be expressed as one boson operator a\, under the condition The 
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transformed form of the Hamiltonian Q can be written as: 



H = SHS- 1 = y2a klk2k ,(at) kl (ai) k *(j--af ai -— a+<j-) ki 
ki ' ' P P 

^E^A^)' 1 ^)^'- f «x«i - + (35) 

U P P 

' ' p p 

mi 

<r-Y,S ni n 2 ,nA4) ni M n2 U ~ 1 *Ui - —ir+<r-) n *. 
^1 • p p 

Note that £3,^3, 7773 and 713 are eliminated by using the conditions given in l|llj) The difference between and 135(1 
is that while in the first the total number of ai and 02 bosons characterize the the system, in the later it is only the 
number of 02 bosons that characterize the system. Therefore the representation is characterized by a fixed number 
j and in ((35(1 . the Hamiltonian is expressed in terms of one boson operator a\. The transformed Hamiltonian H, in 
the Bargmann-Fock space, which play an important role in the quasi-exact solution of the equation Q). It can be 
transformed in the form of the one dimensional differential equations in the Bargmann-Fock space when the boson 
operators are realized as 

ai = — , at = x. (36) 
ax 

The basis function of the primed generators of the system is two component spinor; 

/ x° x 1 ■ ■ ■ x n \ 

Pn,m{x) = I n 1 m ) • (37) 

Y X , X , • ■ ■ , x" 1 ) 

Action of ((35(1 on the ((37(1 . in the Bargmann-Fock space can be written as: 

HP n (x) =J2 P n,m( E )( xn ^ m ) (38) 

The wavefunction is itself the generating function of the energy polynomials. The eigenvalues are then produced by 
the roots of such polynomials. If the E n ,m is a root of the polynomial P n +i, m +i{E) : the series 1(38(1 terminates and 
E n ,m belongs to the spectrum of the corresponding Hamiltonian. The eigenvalues are then obtained by finding the 
roots of such polynomials. 

T-TRANSFORMED HAMILTONIAN 

The constant of motion N characterize the system can be transformed, by the operator T is given by 

N' = TNT^ 1 = (s + pe)a+ ai + pa+a 2 + prja+a^ + r. (39) 

The Hamiltonian (JJJ can be characterized, in the transformed case, e = — s/p and rj = —2r/p. Thus according to 139(1 
the representation is characterized by a fixed number <J2 = j ~ 1- Therefore the transformed Hamiltonian includes 
one boson operator a\, when the condition 1 (11(1 is taken into consideration, it can be written as: 

k ^3 P P 

(roE^AAK)' 1 ^)^' -y- V«i - -*+°-) ei + (40) 
u 3 p p 

<J+ V 7mi m 2 , m4 (a^) mi (ai) m2 (j - — - - -v+°-) mi + 

a — / ' m 3 p p 

mi 

<7-VX n 2 ,nMt) ni (aiT 2 (j - — - f a+ai - -v+o-) n \ 

n 3 P P 
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The Hamiltonian can be expressed as one dimensional differential equation in the Bargmann-Fock space. 

In order to obtain exactly or QES Hamiltonians one can use the same basis given in 1)370 ■ The physical Hamiltonians 
can be obtained and solved by the choice of the appropriate values of the s and p. Consequently we have obtained 
two classes of Hamiltonians whose spectrum can be obtained (quasi)exactly. In the following section we discuss the 
symmetry properties of the general Hamiltonian Q . 

DISCUSSION: THE UNIFIED APPROACH 

In the previous sections transformations of the general two-mode bosonic 2x2 Hamiltonian were treated and its 
differential realization in the Bargmann-Fock space have been obtained. In addition, we note that our approach 
unifies various Lie algebraic methods. The algebraic approach to the finite dimensional part of spectrum consists in 
expressing the Hamiltonian, in terms of the generators of an Lie algebra. When the Hamiltonian can be written in 
terms of the Casimir invariants of the algebraic structure then the eigenvalue problem Hip = Eip can be solved in the 
closed form, giving rise to energy formulas |28t l29j . Otherwise, in general, the spectrum of H cannot be calculates in 
the closed form. If it is written in terms of the bilinear combinations of the generators of the Lie algebra then the 
eigenvalue problem Htf) = Eip must be solved numerically or one can obtain quasi exact solutions. 

A convenient way to construct a spectrum generating algebra is to introduce the boson representations of the 
corresponding Lie algebra. Let us turn our attention to the purely bosonic part of the Hamiltonian Q . This can be 
obtained by setting (3 — 7 = 5 = 0. Under some certain conditions one can obtain and Hamiltonian which include 
bilinear products of the bosonic operators: 

01, a~2 a>2, Oi a>2, at a>i (41) 

generates the Lie algebra su(2) and one can recast these four generators in a more familiar form by introducing the 
three generators 

J = — (afai — a-2 a z) > J+ = a t a 2i J- = a^ai- (42) 

The one variable differential realizations of the generators 1)421 can be obtained by the transformation procedure 
given in the previous sections and they play an important role to the QES of the Hamiltonian which underlying su{2) 
symmetry. 

The other important Lie algebra that included by the Hamiltonian Q is the sw(l, 1) algebra. Bosonic representa- 
tions of the sit(l, 1) algebra is given by 

Kq = — (afai + a\a2 + l) , K + = a^a^ , K- = a 2 a\. (43) 

Therefore by the appropriate choice of the parameter on j \ of one can construct the Hamiltonian possesses the 
symmetry of the su(l, 1). The transformations of the operators of the su(l, 1) lead to the one dimensional differential 
realizations of the algebra. Both in the su(2) and su(l, 1) algebras the transformed operators are in the form of 
the operators of the si-algebra. One can treat bound state problems by using the su{2) algebra and scattering state 
problems su(l, 1) algebra. Thus one can have transition from one bound state to another and from a scattering state 
to another. In order to calculate from bound states to scattering states we need a larger 

algebrajUHEl- We can 

construct this algebra by considering the bilinear combinations of the bosonic operators 

a^ai, a^ei2, CL2, ajai, afoj, 0,20-1, aifii, 02^2, a i a i 1 a 2 a 2- (44) 

One can show that these 10 operators close under the sympletic algebra Sp(4, R). The algebra Sp(i, R) contains both 
bound state algebra su(2) and scattering state algebra su(l, 1). Therefore the algebra provides a unified treatment 
within both bound and scattering states. Its generators can connect all states in the same potential. It thus provide 
a unified approach to the one dimensional problems. 

In addition to the Lie alge bras su(2),su(l,l) and Sp(i,R) the Hamiltonian QJ include two other Lie algebras; 
osp{2, 1) and osp(2, 2) [32, H3- The natural step to relate the Hamiltonian and osp(2, 1) algebra is to express the 
Hamiltonian as linear and/or bilinear combinations of the operators of the osp(2, 1) algebra. The algebra su(2) can 
be extended to the osp(2, 1) with the operators: 

V+ = (T + a 2 , V- = -o+ai, W+ = o-a\, W- = o-a\ (45) 



8 



or by introducing the operators 

V + = cr_a 2 ,K. = -a-di,W + = a+af, W- = o+a\ (46) 

The transformation of the operators of the osp(2, 1) by the operators S or T gives its one variable 2x2 matrix 
realizations which are useful for practical applications. 

One of the major symmetry group candidates for spin one-half particles is the supergroup osp(2, 2) which has four 
even and four odd generators. Its even generators can be represented by bosons while odd generators are represented 
by combinations of the fermions and bosons. The superalgebra osp(2,2) might be constructed by extending su(l, 1) 
algebra with the fermionic generators. It is possible to express two set of fermionic generators to extend the su(l, 1) 
algebra to the osp(2, 2) algebra. These are given by 

V+ = ct-(4, V- = a-ai, W+ = <T+a^ 7 W- = <r + a 2 (47) 
V+ = a + a% ,V- = a + ai,W + = a-a£ ,W- = a-a,2. (48) 

Among these let us also mention here, the procedure given here can also related to the nonlinear Lie algebras that 
have been great deal of interest because of their several significant applications in several branches of physics. Let us 
illustrate these relation on an example. Consider The effective Hamiltonian, which represents the Jaynes-Cummings 
model with Kerr nonlinearity, have been expressed as|34j| 

H = toa + a + — uuq<jq + «(a + c_ + aa + ) + Xa + aa + a (49) 

where n and A are coupling constants of the field and atom and coupling constant of the field and Kerr medium 
respectively. The Hamiltonian H49fl can be expressed in terms of the Hamiltonian Q by an appropriate choice of 
parameters and the number operator of this structure becomes: 

N = a+ ai + icro- (50) 
The invariance algebra of the Hamiltonian Ij49(l is generated by introducing the generators, 

Y + = cii<j + , Y- = a^(T_, Yq — a^ai + <jq (51) 

yields the commutation relation: 

[Y , Y±) = ±Y±, \Y+,Y-] = (1 + 2Y )(Y - N) - i (52) 

The commutation relation [Y + , YL] is a polynomial in Yq. Since the deformation is quadratic in Yq, we have a 
quadratic algebra. The algebras of type i(52"|) have been considered as deformed su(2) algebra. We can easily express 
the Hamiltonian (|49(l in terms of the generators of the deformed su(2) algebra, 

H = uj(2N - Y Q ) + u (Y -N) + k{Y + + F_) + A(2iV - Yq) 2 (53) 

Note that the number operator N is associated with the conserved quantity of the physical system and it commutes 
with the generators Y+, K_ and Yo- Consequently, in this article, we have shown various Lie algebraic approaches can 
be treated in a unified framework. 



CONCLUSION 



In this paper we have prepared a general method to obtain the solution of two boson and one fermion Hamiltonian. 
By using either solution of number operator or similarity transformation, we have been able to provide a QES of the 
various physical Hamiltonians. Furthermore, it has been given that two boson Hamiltonian can be reduced to single 
variable differential equation in the Bargmann-Fock space. It is also important to mention here that the methods 
given here can be used to solve higher order differential equations. 

The algebras can be realized in several ways. For practical applications, the realizations are given in terms of the 
boson creation and annihilation operators. In this paper we have discussed the connection between two-boson and one 
variable differential realizations of the various Lie algebras. The realization in terms of the one variable differential 
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equation directly leads to the usual Schrodinger formulation. The method given here, can easily be extended to solve 
the Hamiltonians that include multi-boson or multi fermion-boson systems. We have presented a first step toward 
the extension of the formulation to obtain solution of the various physical problems. Finally, we have shown that the 
technique given here provides a unified approach to the various algebraic techniques. 
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